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This note describes fitting formulae for the gravitational waveforms generated by a 



^^ ' rapidly rotating neutron star (e.g., newly-formed in the core collapse of a supernova) as it 

ps| ■ evolves from an initial axisymmetric configuration toward a triaxial ellipsoid (Maclaurin 

(n: 

spheroid ^ Dedekind ellipsoid). This evolution is driven by the gravitational radiation 



r^ \ reaction (a special case of the CFS instability). The details and numerical results can 
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Q"^ ! The waveform (including the polarization) is given by Eq. (3.6) of LS. Since the 



found in [Lai & Shapiro, 1995, ApJ, 442, 259; Here referred as LS]. 
I will use the units such that G = c= 1. 



waveform is quasi-periodic, I will give fitting formulae for the wave amplitude h (Eq. [3.7] of 
LS) and the quantity (dN/dlnf) (Eq. [3.8] of LS; related to the frequency sweeping rate), 
from which the waveform hj^{t) and hx{t) can be easily generated in a straightforward 
manner. 

Wave Amplitude: The waveform is parametrized by three numbers: fmax is the 
maximum wave frequency in Hertz, Mi. 4 = M/(1.4M0) is the NS mass in units of 1.4M0, 
Rio = i?/(10km) is the NS radius in units of 10 km. (Of course, the distance D enters the 
expression trivially.) It is convenient to express the dependence of h on t through / (the 



^ This note was written in 1996. It was not intended for publication. Since I have 

been getting requests from people interested in GW data analysis about the waveform 

information, I thought it might be useful to put this note on gr-qc, so that I don't have to 

spend time looking for the TeX file every time I get a request. 
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wave frequency), with f{t) to be determined later. A good fitting formula is 



h[f{tyjma.,M,R] = -—A(-L- 
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where 



4 _ / (/max/1756)2-^, for /^a^; < 500 Hz; ,,Ht}i f =f r3/2,,-i/2 

" I (/^a./1525)3, for /_. > 330 Hz; ^^^^ ^"^''^ " ^"^"^^^lo '^'^-^ ' 



Note that if we want real numbers, we have 
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Number of wave cycles per logarithmic frequency: The fitting formula is 
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Note (i): using equations (l)-(2), we obtain the "characteristic amplitude": 
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which agrees with Eq. (3.12) of LS to within 10% [Note that in Eq. (3.12) of LS, the factor 
/^/^ should be replaced by {R^q M^ J /)^^^, similar to the above expression.] 

Note (ii): The accuracy of these fitting formulae (as compared to the numerical 
results shown in LS) is typically within 10%. When / is very close to fmax, the error in 
the fitting can be as large as 30%. 

Note (iii): The frequency evolution f{t) is obtained by integrating the equation 
f^/ f = —\dN/d\nf\ (note that the frequency sweeps from fmax to zero). For example, we 
can choose t = at / = 0.9 fmax- (Note that one should not choose t = at / = fmax as the 
time would diverge — the actual evolution near fmax depends on the initial perturbations) . 



Once f{t) is obtained, the waveform can be calculated as (cf. Eq. [3.6] of LS): 

\ -L- cos^ 9 
h+ = h[f{tyj^ax,M,R]cos<^{t) , 

hx = h[f{tyjmax, M, R]sm<5{t) cose, 

where 9 is the angle between the rotation axis of the star and the line of sight from the 
earth, and $(t) = 27r J f{t)dt is the phase of the gravitational wave. 

Note (iv): fmax typically ranges from 100 Hz to 1000 Hz (see Fig. 5 of LS); Mi. 4 
and Rio are of order unity for realistic neutron stars. 



